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I. INTRODUCTION 


The purpose of this research is to develop high- and low- frequency meth- 
ods to analyze various radiating elements located on aerospace vehicles with 
combinations of conducting, nonconducting, and energy absorbing surfaces 
and interfaces. Current methods of analysis require limitations on the geo- 
metrical shapes used to model the structure and on the surface conductivity 
of the outer skin. In addition to being computationally inefficient, most ex- 
isting vehicle antenna simulation codes allow surface impedance to be mod- 
elled as perfectly conducting only. While this may be sufficient for many 
configurations, it will not adequately model future airframe designs which 
incorporate advanced composite materials. The performance of on-board 
systems can be significantly affected by the electromagnetic characteristics 
of these materials. If analysis tools to accurately predict the radiation char- 
acteristics of antennas located on these advanced material airframes were 
available, antenna and skin materials could be synergistically designed and 
located to actually enhance performance. However, at present, the analyt- 
ical and simulation tools necessary to take advantage of such designs are 
not available. 

During this report period, the research focuses to develop fundamental 
concepts, techniques, and algorithms which will remove some of the present 
limitations in predicting the radiation characteristics of antennas on com- 
plex aerospace vehicles. In order to accomplish these tasks, the following 
subjects are examined. A first major area of investigation is to develop rig- 
orous analysis of surface discontinuities of metallic and nonmetallic surfaces 
using the equivalent surface impedance concept and Green’s functions. Sec- 
ondly, the effects of anisotropic material on antenna radiation patterns are 
examined by the equivalent surface impedance concept which is incorpo- 
rated into the existing numerical electromagnetics computer codes (NEC[l] 
and ESP [2]). Measurements techniques have also been developed to charac- 
terize the electrical properties at microwave frequencies. Samples of typical 
materials used in airframes have been tested and the results are included 
in this report. Finally, the fundamental concepts of Precipitation Static 
(P-Static), such as formulations and analytical models, are addressed. A 
computer code has been used to model the P-Static process on a simple 
structure. 
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II. COMPOSITE MATERIALS 

A. INTRODUCTION 

In this chapter, analytical methods to predict the effects of composite ma- 
terials on radiation and scattering patterns, and measurement techniques 
to determine the electrical properties (permittivity and permeability) of 
composite materials are discussed. First, the equivalent surface impedance 
concept outlined below was incorporated into the NEC code to obtain 
backscattering radar cross sections of a partially coated PEC cube and 
radiation patterns of a quarter wavelength monopole antenna mounted on 
a finite PEC plate coated with composite material. Second, progress in the 
Green’s function moment method analysis to treat dielectric surface dis- 
continuities is described. Finally, measurement techniques used at ASIJ’s 
EMAC facility to determine electrical properties (i.e., permittivity and per- 
meability) of composite materials are discussed. 

B. EQUIVALENT SURFACE IMPEDANCE 

1. INTRODUCTION 

The use of nonmetallic materials in aerospace vehicles, such as helicopters, 
are increasingly high because they increase the strength, enhance the per- 
formance of antennas, and modify the radar cross section of airframes. 
However, the analysis of the EM characteristics for such complex struc- 
tures are difficult except for certain types of canonical problems. Even 
though there are many electromagnetics computer codes available, most 
of them (such as NEC, ESP, and AARPC[3]) are developed for structures 
consisting of only perfectly conducting material. In this section, the treat- 
ment of the perfectly conducting objects coated with composite material 
using the NEC and the ESP codes is discussed. 

To treat the PEC structures coated with composite material, the Equiv- 
alent Surface Imped since approach is adopted, and then the resulting impedance 
is incorporated into the NEC code. To verify the present approach, some 
numerical examples of antenna and scattering problems are considered, and 
comparisons are made between measured and calculated data. 
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2. THEORY 


In developing the equivalent surface impedance for the PEC plate coated 
with composite material, two assumptions are made. The first assumption 
is that the plate structure is of infinite extent. The second is that the 
plane wave is normally incident upon the interface between the air and the 
composite material. By doing these, we are neglecting the edge effect due 
to the finite size and the effect of the oblique incident angle. A two-sheet 
impedance approximation is briefly described below. 

Let us consider the backscattering of a perfectly conducting finite plate, 
whose top side is coated with a composite material, illuminated by the plane 
wave normally incident on the air-material interface as shown in Figure 1. 
The material has a permittivity of e, a permeability of /i, and a thickness 
of t. Now the original problem shown in Figure 1 is replaced by the two 
electric surface currents, Jt at the top surface and J B at the bottom surface, 
as shown in Figure 2. 

Application of the boundary conditions for the tangential electric field 
at the interfaces shown in Figure 2 (i.e., St and Sb ) yields 

E i B + E t B (J T ) + E t B (J B ) = 0 (1) 

at the bottom surface, and 

E' t + El(J T ) + E't ( Jb) = E[ t (2) 

at the top surface, where subscripts T and B denote the top and bottom 
surfaces, respectively. The superscripts i, s , and t represent the incident, 
scattered, and total fields, respectively. E\ T is the total electric field just 
inside of the air-material interface of the original problem illustrated in Fig- 
ure 1. Using the reflection coefficient for the plane wave normally incident 
upon the air-material interface, the total electric field E\ T can be written 
as 

E\ t = E't + Ej(Jb , Jt) = (1 + (3) 

where the reflection coefficient T is given by 

_ tytanhfrt) -jo ^ 

i i tanh( 7 <) + 1 ] 0 " 
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Figure 1: A geometry for the bacskcattering of perfect electric conductor 
(PEC) coated with a composite material of permittivity t and permeability 
H due to the normal plane wave incidence. 
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Figure 2: Equivalent problem of Figure 1 for the two sheet impedance 
approximation. 
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In (4) 7 is the propagation constant in the composite material. If we use the 
current sheet approximation [4] for fields radiated by the electric currents 
( Jj and Jb) placed on the top and bottom surfaces of Figure 2, the total 
electric field just inside the air-material interface is obtained as 


E\ t = 


(1 + r)sin(kt) 
L**2Mn(fcf)+i(l + r)ei* 


Jt 


( 5 ) 


where k is the free space wave number, and 7 and 7 0 are intrinsic impedances 
of the material and the free space, respectively. The term inside the bracket 
of the above equation can be interpreted as an equivalent surface impedance 
which relates the total tangential electrical field to the tangential electric 
surface current at the top surface. Then the equivalent surface impedance 
(Z) can be written as 


7 - Si _ n (1 r)sin(fct) / 6 x 

Jt 2 sin(fct) + j(l + r)e jfc< 

Note that the impedance sheet should be placed at x = 0, and the PEC 
plate needs to be located at x = — t to model the original problem of Figure 
1 using the NEC. 

In the next section, the two-sheet surface impedance is incorporated 
into the NEC code to calculate antenna and scattering patterns. 


3. NUMERICAL EXAMPLES 

The equivalent surface impedance of the PEC-backed composite material 
of infinite extent, with a thickness of 0.05”, a relative permittivity of t r = 
10.2 — .7 3.8, and a relative permeability of /i r = 2.12 — j'1.5, obtained using 
(6), is 

Z = 37.04 -;96.14ft. (7) 

Now this impedance value is used to calculate some scattering and antenna 
patterns of finite size structures. 

Figure 3 displays the backscattering RCS of a PEC cube whose one side 
is partially coated with a composite material, whose relative permittivity is 
e r = 10.2 — j'3.8, and relative permeability is fi T = 2.12 — jl.5. The pattern 
is obtained using the NEC when the cube is illuminated by a soft polar- 
ized plane wave. The cube has side dimensions of 1.982” and the operating 
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frequency is 10 GHz. The dimensions of the material patch are 0.992 by 
0.992" and the thickness is 0.05". The segment length used for the NEC 
wire- grid model of the partially coated PEC cube is 0.15 A and the total 
number of segments is 1,812. The CPU time on the CRAY X-MP/18se su- 
percomputer to obtain the complete backscattering RCS pattern is 31,741 
seconds. The calculated NEC results are compared in Figure 3 to the mea- 
sured data obtained at ASU’s EMAC facility and the calculated data using 
the FD-TD method. The results agree rather well near the normal direction 
to each face of the cube. However, some discrepancies exist in scattering 
directions away from the normal to each face. There are several possible 
reasons for these discrepancies. In an attempt to resolve the differences, 
the following three items need further investigation. 

• To use a smaller segment length than 0.15A. Presently, we can not 
handle the NEC wire-grid model using a segment size smaller than 
0.15 A because of our limited memory access on the CRAY X-MP/18se 
supercomputer (we are in the process of expanding the memory access 
on CRAY X-MP/18se). 

• To examine the effect of finite size into the equivalent surface impedance. 

• To investigate the angular dependency of the equivalent surface impedance 
(directional impedance concept). 

In Figure 4, the backscattering RCS pattern of a 2A by 2 A PEC-backed 
composite material, which has the same electrical properties as that of 
Figure 3, is presented. As one can see, the material used for the coating 
reduces the RCS at the coated side very significantly (approximately 13dB 
down) compared to the PEC side. From this scattering example, it is 
evident that the composite material is very lossy. 

To see the effects of the same material on the radiation characteristics 
of an antenna, the radiation patterns of a quarter wavelength monopole are 
calculated for three different cases, and they are plotted in Figure 5. They 
are: 

• A/4 monopole on 2 A by 2 A PEC plate. 

• A/4 monopole on 2 A by 2A partially coated PEC plate. 
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• A/4 monopole on 2 A by 2 A fully coated PEC plate. 

The predicted radiation patterns of the quarter wavelength monopole 
for the three different cases do not show significant changes due to the use 
of the composite material. This is in contrast to the scattering patterns 
which were greatly influenced by the presence of the material. There exists 
two possible reasons for minimal changes in radiation patterns for the three 
different cases. The first reason may be that the direct field of the source 
element (quarter wavelength monopole) is dominant in the total pattern 
so that the field reflected from the ground plane does not contribute sig- 
nificantly. The second possible reason is that the radiation from the A/4 
element toward the ground plane is at oblique angles for which our theory is 
not yet very accurate. Experimental verifications of these problems will be 
performed in the next period using measurements at ASU’s EMAC facility. 

4. FUTURE WORK 

As was mentioned in the previous section, the equivalent surface impedance 
given by (6) was derived based on two assumptions. First, we assume 
that the PEC-backed composite material is of infinite extent. The sec- 
ond assumption is that the plane wave is normally incident upon the air- 
material interface. Presently, corrections to both assumptions are under 
investigation. Preliminary results indicate that the new equivalent surface 
impedance will be nonuniform, and its value will be a function of position 
along the material and will also depend upon the angle of incidence. 

Implementation of this new surface impedance approach to the NEC 
and ESP codes will be performed in the next period. Some radiation char- 
acteristics of antennas mounted on complex aerospace vehicles, such as 
helicopters, will also be considered. In addition to numerical predictions 
using the EM codes, measurements of antenna and scattering patterns will 
be made at ASU’s EMAC facility and will be compared with predictions. 
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Figure 3: Backscattering RCS of PEC cube whose one side is partially 
coated with a composite material of e r = 10 — ^3.8 and /i r = 2.12 — j 1.5. 
Pattern cut: x-z plane. Soft polarization. 
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Figure 4: Backscattering RCS of 2A by 2A PEC plate coated with a com- 
posite material of e r = 10 — jZ.S and /i r = 2.12 — j 1.5 using NEC. Soft 
polarization. 
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Figure 5: Radiation patterns of A/4 monopole mounted on PEC plate 
(solid line), partially coated PEC plate (w=lA) (dashed line), and fully 
coated PEC plate (w=2A) (curve formed by black dots) using the NEC. 
c r = 10 — j3.8 and fi r = 2.12 — j'1.5. 
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C. ANALYSIS OF DIELECTRIC DISCONTINUITIES 
UTILIZING GREEN’S FUNCTIONS TECHNIQUES 


1. OBJECTIVE 

The objective of this work is to develop numerical techniques that will 
allow existing codes to analyze dielectric-coated surfaces. Recent trends 
in aerospace vehicle technology point toward extensive use of composite 
materials for major parts of aerospace vehicle structures. Some of these 
composites have anisotropic permittivity as well as anisotropic conductiv- 
ity. In previous reports, we introduced a spectral domain Green’s function 
approach to calculate the currents and the fields on a perfect conductor 
with a dielectric coating of thickness t as shown in Figure 6. Here, we 
extend this work to formulate the surface impedance due to these disconti- 
nuities. This surface impedance for the discontinuity can then be used as a 
correction factor to the surface impedance of an infinite surface to obtain 
the surface impedance of the finite size coating. The impedance of the fi- 
nite size coating is then used to examine the influence of the coating on the 
electromagnetic characteristics of radiation and scattering patterns from 
helicopter structures. We have also explored a spatial domain approach to 
find the surface impedance correction factor. A comparison between the 
two approaches will be described. 

2. EQUIVALENCE THEORY OF DIELECTRIC DISCONTI- 
NUITIES 

There is a considerable literature devoted to the analysis of metallic dis- 
continuities, including codes that are currently used for the analysis of he- 
licopter surfaces. However very little work has been reported for the anal- 
ysis of dielectric discontinuities. Because the discontinuity analysis must 
be valid over a wide range of frequencies, dielectric thicknesses, and dimen- 
sions, the choice of analysis is limited to full-wave techniques. Virtually no 
rigorous full-wave analyses have been reported for dielectric discontinuities. 
Therefore new concepts have to be developed to study these discontinuities. 

Our new concept for treating dielectric discontinuities is shown in Figure 
7, where the finite dielectric coating of Figure 6 is replaced by the difference 
between two coated surfaces. The first surface (surface “a” in Figure 7) is an 
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infinite coated surface with a well known surface impedance. The second 
surface (surface w b”) is a partially coated surface, whose analysis is the 
subject of our investigation to find its equivalent surface impedance. The 
equivalent surface impedance of Figure 7(b) will be used as a correction 
factor to the surface impedance of the infinite structure to find the surface 
impedance of the partially coated PEC. This equivalent impedance will be 
found by studying the boundary conditions at a dielectric-air discontinuity 
in surface w b’\ 

Assuming no charge on the surface and that the dielectric material 
is magnetically isotropic, all fields will be continuous except the normal 
electric field components £ n i, E n 2 (see Figure 8). Using equivalence theory, 
this discontinuity can be represented by an equivalent source of electric 
field equal to the difference in normal field components. Assuming that the 
discontinuity exists at x = x 0 as shown in Figure 8, then the discontinuity 
source can be written as 

E. = ♦ «(* - *.) (1) 

trx ^rx 

where e TX is the dielectric constant of the dielectric material in the x direc- 
tion, * denotes the convolution of two functions. Using convolution theory, 
(1) can be written in the spectral domain as 

E t = (2) 

^rx 

where * denotes the spectral domain with respect to the coordinates x, y. 
The normal field component E n is related to the ground plane current J a 

by 

E nl = x -T t (z)T g (3) 

where A denotes the unit vector, z is height of the observation point, and Zt 
is the transmission submatrix (see the second quarterly report). From (2) 
and (3) 

E e (z) = • t t (k x , k v , z)T g (k x , k v ). (4) 

trx 

The surface fields and currents E,, J t due to this discontinuity equivalent 
field can be found as 

E, = J'T)(k x ,k v ,t- z)E e (z)dz, 
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J, = J^7j(k x ,k v ,t - z)E e (z)dz. (5) 

Finally the correction factor Z, due to dielectric discontinuity can be found 
as . 

Z, = E.Z~ • (6) 

To implement this formulation into the NEC code the currents on the 
ground plane are calculated by first assuming the dielectric coating is of 
infinite size. This current needs to be transformed into the spectral domain 
using one of the commonly used FFT subroutines. Then the discontinuity 
surface impedance is evaluated and transformed back to the spatial domain 
using IFFT. Finally the new surface impedance can be calculated as a 
combination of infinite and discontinuity surface impedances. We hope 
that only one additional rerun of the code will be required. However an 
iterative technique can find an accurate solution. 

As mentioned previously, we also considered an alternative spatial do- 
main approach. The main attractive feature of the spatial domain technique 
is its direct compatibility with the available computer code. However, the 
spatial domain approaches suffer from the following disadvantages: 


1. The Green’s functions corresponding to the transmission functions 
Z, Y, T are very complicated and the extension to anisotropic materi- 
als is not promising. For example, the normal electric field E n at the 
discontinuity can be expressed as 


<■- II 


^ W ^ (1 + _ 1 


2irr 

jk sin 2 (0) 

2 

jk cos 2 (9) 

2 


) 


✓ 1 1 .sin 2 ^)., 1 * 

^ + Jlr “ I»r» ^ ~ ^ r * + 7 Or* 

(1 + 7 ^- 7 ^))] 


j kr k 2 r 2 

where 

6 = tan -1 (- — — ),a = tan _1 (— ). 

y-y o r 

The subscript o denotes the observation point and 


r = y/(x - x 0 ) 2 - (y - y 0 ) 2 + (z - z 0 ) 2 . 


( 7 ) 

( 8 ) 

( 9 ) 
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2. The spatial domain approach is directly compatible with the codes 
but will probably be computationally slow due to the larger number 
of required integrations. 

3. The spectral domain approach can handle more than one disconti- 
nuity. We will at least need to calculate the effects of four different 
discontinuities. 

3. FUTURE WORK 

In the near future we will concentrate on spectral domain approaches to 
modify the computer code and compare predicted results with experimental 
results to verify our new concepts. The modified code will also be used to 
study the radiation and scattering from complex aerospace vehicles, such 
as helicopters, including dielectric-coated sheets. If successful, we will pro- 
ceed to study anisotropic and composite materials using the same spectral 
domain approach. The spatial domain method may be also considered, 
especially in the simple case of isotropic materials. 
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Figure 6: Finite dielectric- coated PEC. 
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Figure 7: Equivalent surface of a finite dielectric coated surface. 
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D. MEASUREMENT OF ELECTRICAL PROPERTY 
OF MATERIALS 

1. INTRODUCTION 

There are many design and analysis problems in the field of telecommuni- 
cations in which it is vital that the precise electrical properties of materials 
be known. The characteristics of a material that are of most importance to 
the electrical engineer are the conductivity (<x), the complex permittivity 
(e r = e' r —je"), and the complex permeability (/i r = #4 — JV-")- The permit- 
tivity of dielectrics is a pivotal quantity in the design of microstrip circuits 
and devices, resonators, capacitors, rf cables, etc.. In some applications, 
the permittivity and the permeability are of great importance, as in the 
design of radar absorbing materials. And with the increasing utilization 
of light-weight advanced composites and low radar-observable materials, 
simple conduction modeling is inadequate to accurately predict the electro- 
magnetic scattering of modern aircraft or the radiation characteristics of 
associated antenna systems. 

The first permittivity and permeability measurements performed at 
ASU’s ElectroMagnetic Anechoic Chamber (EMAC) facility were conducted 
using the Transmission/Reflection (TR) technique. The basic theory and 
the particular implementation of the technique will be discussed. Some 
measurement results will be presented. A number of possible future im- 
provements and goals will also be discussed. 

2. THEORY 

The following explicit equations for the complex permittivity and perme- 
ability from measured broad-band transmission and reflection scattering 
parameters were originally derived by Nicolson and Ross [5] and later by 
Weir [6]. Other formulations have been developed more recently [7], [8], but 
not with the simplicity and elegance of the original. 

Consider an infinite dielectric slab of thickness d as shown in Figure 9. 
Assume the following conditions: 

1. The front and back surfaces of the slab are flat and parallel to each 
other. 
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2. The slab is homogeneous, isotropic, and has permittivity e = Crto and 
permeability fi = fi r n <>, where e 0 and are free space values. 

3. A plane wave is normally incident on the slab at z=0. 

The system shown in Figure 9 can be considered conceptually as three 
cascaded two-port networks. The first two-port consisting of just the air- 
to-dielectric interface, the second consisting of just a section of dielectric of 
length d, and the third consisting of just the dielectric-to-air interface. The 
first and third two-port networks have scattering parameters, by definition, 
in terms of only the reflection coefficient T of a single interface, given by 

r - Izlli - y*_ A . (i) 

r ? + T 7‘ y/^ + 1 

where tj is the impedance of the dielectric and i / 0 is the impedance of free 
space. The second two-port network has scattering parameters involving 
only unidirectional propagation through the dielectric, given by 

T = t~^ d = e -J «v"‘ r<r<i (2) 

where 7 is the propagation constant of the material. By performing the 
appropriate algebraic operations on the three individual two-port scattering 
matrices, it can be shown that the scattering matrix of the system shown 


in Figure 9 is 

r (i-r 2 )r (i-r J )T l 

1 -r 

[$(«)] - 

( 3 ) 

<*-?) I (*-££ 

Solving for T and T in terms of the measured scattering parameters Sn(u;) 

and S 2 iO*>) yields 

r = x±Vx^T 

( 4 ) 

where 

. . 

( 5 ) 

and 

25n(a>) 

Sn(ui) + Sai(u>) — T 
_ l-|5„(u.) + 5 J1 (u.)]r- 

(6) 
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The appropriate sign in (4) is chosen such that |r| < 1 . If we define from 

( 1 ) , „ 

ii+n _p,_ m (7) 


and from (2) 


then 


-LSi-®]*-— 


( 8 ) 


( 9 ) 

( 10 ) 


fl r — y/ClC] 

e ' = vf 

Thus the complex permittivity and permeability can be obtained from the 
measured Sii(w) and S 2 i(w) of a material slab. 


3. IMPLEMENTATION 

The conditions assumed for the TR technique can be achieved by at least 
three physical implementations: free- space method, coaxial air line method, 
and rectangular waveguide method. In the free-space method [9], a large 
flat sample is placed between two spot-focusing horn lens antennas. The 
transverse dimension of the sample must be at least three times the beamwidth 
of the antennas at the focus to reduce the diffracted fields to an insignificant 
level. The free-space method is nondestructive. Because it is not necessary 
to machine the material to fill a sample holder, there are no errors associ- 
ated with air gaps. Errors due to the excitation of higher order modes at 
the air-dielectric interface are also nonexistent. However, the bandwidth 
of a single measurement is limited by that of the antennas, and the free- 
space method requires the fabrication of a relatively large precision (and 
expensive) test apparatus. 

In the coaxial air line method, a sample of material is machined to com- 
pletely fill the cross section of a coaxial air-dielectric transmission line. This 
method yields extremely large bandwidth per single measurement. The up- 
per frequency limit is related to single mode transmission by the dimensions 
of the line. This is a destructive method: a sample of the material must 
be cut and is subsequently useless for its intended application. There is a 
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possibility for error due to air gaps between the Bample and the walls of 
the transmission line. Also, the sample must be cut to an annular shape, 
requiring precision machining operations. 

The waveguide method is very similar to the coaxial airline method. The 
sample is cut to fill the cross section of a rectangular waveguide; however 
a good fit between the sample and the waveguide holder is much easier to 
achieve due to its simple shape. The bandwidth of a single measurement 
is limited to the frequencies between the cutoff of the TEio mode and the 
cutoff of the next higher mode. 

The waveguide method is the implementation utilized at EMAC. At the 
time of this report, we are configured for material property measurements 
in X-band (8.2-12.4 GHz); however the components necessary to extend 
our frequency range from 2.6 GHz to 18.0 GHz are commercially available 
in half-octave bands. 

The S-parameters of the material sample are measured with am HP8510B 
Network Analyzer utilizing the 8515A S-Parameter Test Set and 8341B Syn- 
thesized Sweeper. Phase-stable coaxial cables connect the test ports of the 
network analyzer to coax-to- waveguide adapters. Five inch long waveguide 
sections isolate the material sample from higher order modes excited by 
the adapters. The sample holder is a 0.38 w long 1:1.00 VSWR precision 
waveguide standard (Figure 10.). 

The procedure begins with a full 2-port calibration at the ends of the 
5" waveguides. A sample of material is placed in the holder, the holder 
is inserted between the waveguides, and the amplitude and phase of the 
S-parameters are measured. Because the sample length d is generally less 
than the length of the sample holder, the calibration planes must be math- 
ematically shifted to coincide with the surfaces of the sample. Notice from 
(3) that the scattering matrix (assuming an isotropic material) is symmet- 
ric and reciprocal. By comparing Sn to S 22 s^d S 12 to S 21 , after performing 
the phase shift, a verification of the sample location is possible. Finally, 
the relative permittivity and permeability are calculated from the phase 
corrected Sn and 52 1 (or 522 a^d 5i2). 

There are a number of limitations and drawbacks inherent to this method. 
When the loss tan£ of the material is below about 0.01, the uncertainty of 
ejf and n" becomes unacceptably high, although the method still gives valid 
cj. and fi' r . The equations are unstable for sample lengths of (2n + l)A fl /2, 
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and the electrical length of the sample should be greater than about 10 
to 20 degrees for accurate results. Materials that have a rough surface or 
inhomogeneities can excite higher order modes which lead to inaccuracies. 

4. MEASUREMENTS 

One of the first materials measured at EMAC was Teflon. Teflon is useful 
as a benchmark because its characteristics are well known. The measured 
values of cj. in Figure 11a range from 2.055 to 2.034. Although this par- 
ticular sample was not precisely formed, the measured results are in good 
agreement with values quoted in [8] of 2.044 to 2.048 across X-band. The 
measured relative permeability and /x" are reasonably close to the ex- 
pected values of 1.0 and 0, respectively. 

Measured /x r and e r of a precisely formed sample of nylon 66 are shown 
in Figure 12. The real part of the permittivity is within 2% of values quoted 
in [10]. The measured fi' T is nearly 1.0 as expected for this nonmagnetic 
material. 

The data in figures 13 and 14 represent values of magnetically-loaded 
elastomeric surface wave absorbing materials. No values for comparison 
were immediately available but a n' T greater than unity was measured, as 
expected for a magnetically-loaded material. 

The measured e T and /i r of a 0.025” thick sample of woven Kevlar are 
shown in Figure 15. The uncertainty of this measurement is believed to be 
high due to very thin sample and its fibrous surfaces. The measurement 
still gives an indication of the constituent parameters of the material. The 
fiberglass of Figure 16 is the type used in the construction of modern air- 
craft. Twelve plys of 9 mil thick woven fiberglass layed up in epoxy formed 
this 0.114" thick material. 

5. CONCLUSION 

The Transmission /Reflection technique of measuring the electrical prop- 
erties of solids simultaneously provides complex relative permittivity and 
permeability information using measurement fixtures and mathematical al- 
gorithms that are simple and inexpensive. 

Future goals include improved measurements and variations of the tech- 
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clique and algorithms. A precisely formed sample of Teflon has been con- 
structed (but not yet measured) for future benchmark measurements. Time- 
domain gating will be utilized to greatly reduce the effects of source and 
load mismatch that are residual to the calibration. This will virtually elim- 
inate the ripple that is evident in the results presented here. The TRL 
technique [11] for calibrating the network analyzer will be considered. This 
would allow a much faster calibration using less elaborate standards, while 
maintaining the accuracy of a conventional 2-port calibration. The calibra- 
tion currently takes approximately 90 minutes and degrades significantly 
in a few hours. Alternate measurement techniques [12] utilizing the same 
waveguide fixtures will be considered to provide in-house intermeasurement 
comparisons. Alternative algorithms [8] may provide results when the algo- 
rithms presented here fail (when attempting, for example, to measure very 
high-loss, nonmagnetic material). 
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Figure 9: An infinite dielectric slab of length d. 
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Figure 10: A schematic diagram of the measurement fixture at ASU’s 
EMAC. 
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Figure 12: Measured e r (a) and /j r (b) of nylon 66. 
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Figure 16: Measured c r (a) and n T (b) of fiberglass. 12plys, each ply 9mils 
thick. Sample thickness =0.114". 
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III. PRECIPITATION-STATIC (P-STATIC) 

A. INTRODUCTION 

The term precipitation static refers to a general class of interference prob- 
lems in the communication system of an aircraft, flying through a precip- 
itated environment [13]. The interference in the receiving or transmitting 
signals is due to the ice crystals in the precipitation. The degree of this 
interference can be high enough to block the communication for a long 
time. Studies, mostly experimental, have been undertaken to determine 
the cause of the interference and how it can be reduced to a tolerable level, 
if not completely eliminated. 

From early studies, it was determined that one type of precipitation 
static is caused by corona discharge from the aircraft which is charged by 
a triboelectric charging mechanism. 

Triboelectric charging mechanism occurs when two dissimilar materi- 
als exchange opposite charges while they collide and then separate. One 
material acquires a negative charge, leaving the second material with a 
positive charge. The amount of charge depends on the materials and their 
momentum with respect to one another. In the case of an aircraft flying 
through the precipitated environment with ice crystals, the aircraft is nega- 
tively charged, leaving the ice crystals with a positive charge. The charging 
process on the aircarft, flying through precipitation, can be viewed as a con- 
stant current source which increases the potential of the aircraft. Rise in the 
potential continues until it exceeds the corona threshold potential, causing 
the corona discharge. An equilibrium condition is reached when the total 
corona current becomes equal to the source current. 

Another type of charging mechanism is due to the gas exhaust from the 
aircraft’s jet engine. This is observed when a turbojet aircraft flies in low 
altitude. 

In general corona discharge consists of a series of extremely short pulses 
that contains RF components over a broad spectrum, interfering with the 
communication system. 
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B. CORONA DISCHARGE MECHANISM 

Corona discharge occurs when electrons acquire enough energy between 
collisions to ionize other neutral atoms, thus producing additional electrons. 
Existence of the initial electrons is maintained by the cosmic rays and other 
high energy radiation. The energy, that these electrons need for the corona 
discharge, is provided by a sufficiently high electric field in the vicinity 
of a charged body, such as the body of an aircraft[13]. The action starts 
when the initial electrons move away from the high electric field region ( 
sharp points on the negatively charged aircraft ), as shown in Figure 1, and 
collide with air molecules, which some become ionized, and they produce 
additional electrons and positive ions that are accelerated by the field. This 
avalanche process continues and grows until the electrons reach a region far 
enough from the high field point where the field becomes too low to allow 
ionization by collision, and the electrons are slowed down. 

The slowed electrons attach to other molecules, like oxygen molecules, 
producing 0 2 ions, as displayed in Figure 2. Since these negative ions 
are less mobile than the electrons, they remain stationary during the dis- 
charge process. The negative charged cloud made of these negative ions 
produce an opposing electric field which reduces the original high electric 
field, hence reducing the distance for the next avalanche process. Mean- 
while, the cloud of positive ions left from the first avalanche tend to increase 
the field between itself and the high electric field point, thereby allowing the 
avalanche to be initiated more readily between the cloud and the high field 
point. The next supply of the electrons for the avalanche are provided by 
the photoemission from the high negative points on the aircraft. The dis- 
charge process continues as a series of avalanches with successively shorter 
distances, while the negative oxygen ions cloud approaches the high field 
discharge point. The positive ions cloud are drawn to the high negative 
points and are neutralized. The negative space charge region of the 0 2 
reduces the field at the high negative point (see Figure 3), thus blocking 
further avalanches and discharges until they are gradually swept away by 
the wind and the applied electric field. Then the field on the body rises, 
and a whole new process is repeated. The entire single process is completed 
in 0.2 microsecond at the sea level pressure. 

A carefully prepared laboratory experiment on a negative point corona 


34 


discharge has revealed that the rise times of the pulses are very short, 
less than 0.01 fi seconds at the sea level pressure, and that pulses decay 
exponentially. It has also been shown that both the rise time and the fall 
time are inversely proportional to the pressure, and spatial extent of the 
corona is dependent on the size of the point body. 

C. PROPAGATION OF CORONA DISCHARGE IN- 
TERFERENCE 

Two of the most important aspects in the investigation of precipitation 
static are: 

1. Prediction of points on the body of an aircraft where the corona dis- 
charge are more prevalent. These point are characterized by high 
amplitudes of the normal components of the electric field on the sur- 
face of the aircraft. The relative amplitudes of these fields set the 
intensity of the interference. 

2. The coupling between the various interference points and the antenna 
on the aircraft in the presence of space charge region. In what follows 
attempts are made to investigate the second aspect. 

In order to examine the frequency spectrum of the interference and its 
coupling to an antenna, it is necessary to solve for the transient electric and 
magnetic fields generated by the corona discharge process. Furthermore, 
to closely model the process, these fields have to be solved in an ionized 
air surrounding the discharge points. To implement the latter portion of 
the analysis, one has to use air chemistry formulation to predict the air 
conductivity by ionizing radiation [14]. The ionizing process is modeled by 
the following differential equations: 

^ - 4- V • (n e v e ) = Q + Gn e - an e - /?n e n+ (1) 

^ 1 + V ■ (n+u + ) = Q + Gn e - /9n e n + - 7 n + n_ (2) 

at 

+ V • (n_u_) = Q + a e n e - 7n + n_ (3) 

at 
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The pairs (n e ,t> e ), (n+,r+), and (n_,tL) are the number density and 
velocity of electrons, positive ions, and negative ions, respectively. The 
parameters G, /?, 7 , and a e are the avalanche rate (.s -1 ), the electron- 
ion recombination rate (m 3 /s), and the electron attachment rate (s 1 ), 
respectively. Q(t ) is the ionizing source function due to cosmic radiation, 
and approximately it equals to 10 7 electron-ion pairs /(m 3 s). The velocities 
axe proportional to the electric field, and they are given by 


— # 

V e = -p e E 

(4) 

v+ = P+E 

(5) 

tL = p-E 

(6) 


where p e ,/i+, and //_ are the mobilities of the electron, positive ion, and 
negative ion. G, o e , and p, are functions of the electric field. 

One needs to include the continuity equation to take into account the 
space charge effect. This is given by 



(7) 

Jt = + Jc 

(8) 

Pt — Ps + pc 

(9) 

q e (n + v + — ri-v _ - n e v e ) 

(10) 

Pc — c n — -f 

(11) 


The subscripts s, c, and t stand for ’external source’, ’corona source’, and 
’total source’, respectively. J t and p t become the primary sources for the 
electromagnetic interference. 

The electric and magnetic fields of this interference are obtained via 
solution to Maxwell’s equation: 


Vx£ = -f (12) 

_ - _ dD . 

V x H = J t + (13) 
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V-D = p t (14) 

V-£ = 0 (15) 

Boundary conditions on the surface of the aircraft are required to relate 
the charge and current distribution to the electric and magnetic fields. 

For simplicity, we assume the boundaries are perfect conductors; then the 
relations are 

J, = nxH (16) 

p, = n ■ E (17) 

Equations (1),(2),(3),(7), and (12)-(15) must be solved simultaneously 
to obtain the electromagnetic fields generated by the corona and the tribo- 
electric currents. 

To solve for the transient electric and magnetic fields in the presence 
of space charge region, a finite-difference time domain solution is imple- 
mented. For simplicity, it is assumed that a conducting cube is charged, as 
depicted in Figure 4, and then allowed to discharge through the corona dis- 
charge process. The selection of the simple geometry permits complicated 
computations between the charged bodies and the space charge region to 
be performed. 

Next, the far field patterns of this transient response are analyzed for 
two cases: 

1. In the absence of the space-charged region. 

2. In the presence of the space-charged region. 

Frequency transforms of the far field pattern reveal the spectrum of the 
interference. 

D. TRANSIENT ANALYSIS USING FINITE-DIFFERENCE 
TIME-DOMAIN METHOD 

Finite-Difference Time-Domain (FD-TD)[15] is one of the most suitable 
techniques to analyze the transient electric and magnetic field responses of 
the corona discharge of complex geometries, such as helicopters and aircraft. 

Since the corona discharge behaves like an impulse current source in time, 
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and perhaps in space, its excitation emanates waves over a broad frequency 
spectrum. That sets a limit on the range of frequencies for a given FD-TD 
code. This is mainly due to preset cell size and time- step set used in the 
finite-difference method. 

The frequency spectrum analysis of the excited and transmitted waves 
are also limited by the time-step or the time-sampling frequency. This 
sampling frequency, which is generally referred to as the Nyquist rate, must 
be twice the highest frequency of the waveform. In the present analysis, 
the sampling rate equals to Ar 1 where At = 1.67 x 10 -10 seconds is the 
time-step used in the FD-TD code. Hence, one can determine the frequency 
spectrum of waves up to / 0 = 3 x 10 9 Hertz. 

Figure 5 displays the excitation and the transient response of the cube. 
Also shown are the electric field components, E x , E y , and E z , at the front 
face (YZ- plane) of the FD-TD volume, in the absence of the space-charged 
layer. Normal components of the electric field on the surface of the cube 
were used as sources for the charging process. From the results of Figure 
5, it is evident that the normal component of the electric field ( E x ) has 
different transient response from the transverse components of the electric 
field (E y and E z ). This indicates that the normal component has different 
frequency spectrum than either of the transverse components. Figure 6 
shows the frequency spectrums of the excitation and the electric fields whose 
transient responses are shown in Figure 5. The frequency spectrums of 
all transient responses are well below the maximum allowable frequency 
f 0 = 3 x 10 9 Hz. The differences between the transient responses are mainly 
due to the structure of the cube. 

To see how the space-charge region alters the transient responses, and 
hence changes the frequency spectrums , a model of the space-charged 
layer on the surface of the cube was included. Figure 7 shows the transient 
response of the fields in the Y Z-plane. The space charge layer behaves as 
variable air conductance, and it causes the wave to decay. This can be 
verified by the Fourier transform of the fields in the Y Z-plane. Figure 8 
displays the frequency spectrums of the transient waveforms of Figure 7. It 
is apparent that the spectrums of the electric field components in the space- 
charged medium are different from those of the nonspace-charged medium 
shown in Figure 6. 
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E. FUTURE WORK 

Results based on the simple square model of an aircraft clearly shows how 
the corona discharge can cause spurious signals to propagate. In future 
work, transient response of an aerospace vehicle model, such as helicopter, 
will be examined, so that the wave solution for the aerospace vehicle and 
its resonances can be included. 
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Figure 1: Positions of the neutral molecules and the electrons in the high 
amplitude electric field before the avalanche. 
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Figure 2: .Positions of the positive ions and the negative ions with respect 
to the high amplitude electric field after the avalanche. 
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Figure 3: Positions of the negative oxygen ions after the avalanches are 
completed. 
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Figure 6: The frequency spectrums of the excitation said the electric fields 
at the Y Z-plane, in the absence of the space-charged layer. 





Figure 7: The excitation and the transient response of the cube at the 
KZ-plane, in the presence of the space-charged layer. 
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Figure 8: The frequency spectrums of the excitation and the electric fields 
at the YZ- plane, in the presence of the space-charged layer. 
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